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Abstract. We provide a detailed description of solutions of Curve Shortening 
in R" that are invariant under some one-parameter symmetry group of the 
equation, paying particular attention to geometric properties of the curves, 
and the asymptotic properties of their ends. We find generalized helices, and a 
connection with curve shortening on the unit sphere § n_1 . Expanding rotating 
solitons turn out to be asymptotic to generalized logarithmic spirals. In terms 
of asymptotic properties of their ends the rotating shrinking solitons are most 
complicated. We find that almost all of these solitons are asymptotic to circles. 

Many of the curve shortening solitons we discuss here are either space 
curves, or evolving space curves. In addition to the figures in this paper, we 
have prepared a number of animations of the solitons, which can be viewed at 
http : //www.youtube . com/user /solitons2012/videos?view=l. 



A family c : (to, t\) x E — > W 1 of smooth parametrized curves in R" evolves by 
Curve Shortening if it satisfies 



where s is arc length along the curve, and p a denotes orthogonal projection along 
the vector a, so that p Cs c t is the component of the velocity vector c t normal to 
the curve. If the curves are parametrized so that c t _L c s at all times, then (1) is 
equivalent with 



Curve Shortening for compact curves has been extensively studied, especially in 
the case of curves in the plane or on a surface (see e.g. [6].) Here we look at Curve 
Shortening for noncompact curves in K™. More specifically, we consider solitons, 
i.e. self-similar solutions. Since noncompact curves can be geometrically much more 
complex than compact curves, it is to be expected that Curve Shortening for non- 
compact curves will present a much wider range of phenomena than for compact 
curves. 

In geometric analysis, solitons are considered to be solutions of some geometric 
evolution equation that are invariant under some subgroup of the full symmetry 
group of the equation. The equation (1) is invariant under Euclidean motions in 
space, translations in time, and parabolic dilations of space-time. If 25 is the group 
generated by these transformations of space-time, then we will classify and describe 
all solutions of (1) that are invariant under some one-parameter subgroup of 25. In 
§2 we classify the one-parameter subgroups of 25 up to conjugation and find three 
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essentially different kinds of one-parameter subgroup and thus three types of self 
similar solution of (1) (see Table 1). 

Solutions with static symmetry: The first kind of solution we encounter 
("category A" in Table 1) consist of curves that at each moment in time 
are invariant under some family of Euclidean motions. The most general 
static symmetry that is preserved by Curve Shortening is of the form 

(2) e eM C + 0v, (S6l) 

for any fixed skew symmetric matrix M and any fixed vector v with Mv = 0. 
Rotating and translating solitons: The solutions in category B (Table 1) 
are defined for all time iel. They are characterized by the fact that the 
curve c(t, •) at any time t £ R is obtained from the initial curve c(0, ■) by 
a combination of translation and rotation, i.e. 

Here v is the translation velocity and A = —A 1 is an antisymmetric matrix 
representing the rate at which the soliton rotates. 
Rotating and dilating solitons: Solitons in category C (Table 1) are char- 
acterized by the condition that the curve c(t, •) at any time t is obtained 
from the curve c(0) at time t = by a combination of dilation and rotation, 
i.e. 

c(U) = (l + 2at)^ A c((U), 
where a ^ is a constant and A = —A* is the skew symmetric matrix that 
describes the rotation. For a = we have 

c(t,t) = e tA c(0,0 

instead. 

Some special well-known solutions of Curve Shortening include straight lines 
(which fall under all three categories); circles (which have static symmetry, and 
which also evolve by dilating), the Grim Reaper, i.e. the graph of y = — log cos a; 
(which moves by steady translation upward along the y-axis) ; the Yin- Yang curve 
(a spiral in the plane which moves by steady rotation); the Abresch-Langer curves 
(plane curves that evolve by shrinking dilation 1 ); the Brakke wedge (a convex plane 
curve that evolves by expanding dilation.) See Figure 1. The case of plane solitons 
has been fully analyzed by Halldorsson in [8] , where one can also find many figures 
of plane solitons. 

Summary of main results. The existence of special solutions with static sym- 
metry was first reported in [3, 4, 2] where it was noted that under Curve Shortening 
a helix in R 3 will remain a helix, and that this helix will shrink to its central axis. 
Since Curve Shortening preserves the symmetries (2), any solution of (1) whose 
initial curve is invariant under (2) will at all times be invariant under (2). The 
symmetries in (2) form a one parameter group so that curves invariant under this 
group are exactly the orbits of the group action on R ra . The Curve Shortening flow 
on curves invariant under (2) can thus be reduced to a finite system of ODEs. In §3 



Abrcsch and Langer derived and analyzed the ODE that describes planar dilating solitons. 
A main point of their paper [1] was to show that while most of these curves are not compact, 
there is a discrete family of compact immersed curves among them. Nonetheless, we will call both 
compact and non compact purely shrinking solitons "Abresch-Langer curves." 
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we derive this system and find that after a time change it is equivalent with the first 
order linear system x' = M 2 x, which is trivially solved. Translating the solutions 
of x' = M 2 x back to curves evolving by (1), we find two somewhat different cases. 
For v^O the solutions behave as in the Altschuler-Grayson helix example, namely, 
in backward time they are helix-like curves wound on a cylinder of radius 0(v— t) 
while in forward time the curves converge at an exponential rate to the line through 
the origin in the direction of v. For v = 0, solutions only exist for t € (— oo,T) for 
some T. At time t the curve is contained in a sphere of radius \/2(T — t). After 
rescaling to the unit sphere, the solutions converge both as t — > — oo and as t — > T 
to different great circles. In fact, after changing time, the rescaled solutions are 
solutions to Curve Shortening for immersed curves on the unit sphere S n . 

Starting in §4 we consider solitons. Any rotating-translating or rotating-dilating 
soliton c(t) is completely determined by its initial curve c(0). A curve C is the 
initial curve of a rotating-dilating or a rotating-translating soliton exactly when it 
is the solution of the second order ODE 

(3) C ss = (a + A)C + v + AC s , 

where A is determined by the condition C s _L C ss . The parameters a, A and v are 
as above. 

The case of non-rotating solitons, where A = 0, turns out not to present any new 
curves in higher dimensions. On one hand it is a well-known folk theorem that all 
purely translating solitons are planar, so that up to rotation and scaling they are 
equivalent to the Grim Reaper (we prove a more general statement in Lemmas 5.1, 
5.2, and 5.4). On the other hand, it is also known that all purely dilating solitons are 
planar. They are therefore copies of the Abresch-Langer solitons. Again, Lemma 
5.1 presents a generalization of this statement. 

The classification in §2 of one-parameter subgroups of the complete group © of 
symmetries of Curve Shortening shows that we do not have to consider solitons 
that simultaneously dilate and translate. While the soliton equation (3) certainly 
has solutions when a ^ and v =^ 0, these curves turn out to be non-translating 
solitons whose center of dilation and rotation is not at the origin. We therefore 
only consider rotating-translating or rotating-dilating solitons. 

In studying the global behaviour of the curve we distinguish between unbounded 
ends of a soliton, and ends that stay within a bounded region. We show in §6 
that the unbounded ends can be decomposed into high curvature and low curvature 
parts. The high curvature parts are characterized by the condition that the tangent 
and the vector a = (a + A)C + v are not almost parallel. These high curvature 
arcs are short (their length is 0(||a|| -1 )), and the further away they are from the 
origin, the better they are approximated by a small copy of the Grim Reaper. (See 
Figure 4.) On the remaining low curvature arcs, which comprise most of the soliton, 
the tangent C s is almost parallel to a, so that these arcs are solutions of 

(4) A^^ = (a + .A)C + v+ a small error term. 

as 

After reparametrization this equation is linear, and easily solved. It should provide 
a good description of "medium length" sections of solitons. This description can in 
certain cases be shown to be globally correct (see below, and Lemma 7.7.) 

Of all the rotating-dilating solitons, the ones that rotate and expand are simplest 
to describe. The distance to the origin from a point on a rotating-expanding soliton 
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is a function on the curve that has exactly one global minimum, and no other critical 
points. On both ends of the curve the distance to the origin grows without bound. 
Both ends are low curvature arcs, and using the approximate equation 4 we show 
that the ends are asymptotic to generalized logarithmic spirals (Lemma 7.7). In 
the very special case where .4 = these spirals reduce to straight lines, and the 
solitons are just the Brakke wedges. 

For purely rotating solitons the approximate equation (4) (with a — 0, v = 0) 
suggests that the ends will merely rotate. However we find in §6 that the ends 
must be unbounded, which shows that the ends of purely rotating solitons slowly 
"spiral off to infinity." If the null space N(A) is non-empty then either the soliton 
is contained in a proper subspace of K", or else N(A) contains a line £ such that 
the orthogonal projection of the soliton onto I is injective (Lemma 5.4). This must 
always be the case in odd dimensions; e.g. in M 3 the null space N(A) is the axis of 
the rotation generated by A, and Lemma 5.4 says that the soliton is either planar, 
or a graph over the axis of rotation. 

A similar result holds for translating-rotating solitons (a = 0, Av = 0, 0.) 
Here the projection of the soliton onto its translation axis (the line through the 




Figure 1 . For three dimensional solitons the soliton equation (3) 
contains three essential parameters: the dilation rate a, the rota- 
tion rate u>, and the translation velocity v z . For each fixed choice 
of (ol,uj,v z ) many solitons exist, namely, one through each point 
C G K 3 and unit tangent vector T. This figure shows a few solitons 
in R 3 . The named solitons are well known planar examples. 
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origin in the direction of v) is bijective, or else the function x t— > (v, x) has a unique 
minimum on the soliton. In the latter case the projection onto the translation axis is 
two-to-one except at the point where (v, x) is minimal (see Figure 2.) The quantity 
(v,x) is unbounded along each end. If (v,x) — > oo along an end, then along this 
end the distance to the translation axis must also become unbounded. Any end 
along which (v,x) —¥ — oo converges to the translation axis (Lemma 6.4). 

The most complicated solitons are the rotating shrinking solitons (a < 0, A ^ 0, 
v = 0). First, given any A any circle of radius a that is invariant under 

the rotations e SA is trivially a soliton that dilates and rotates with parameters 
(a, A). Similarly, any Abresch-Langer curve in N(A) is a dilating-rotating soliton 
with the same parameters. In §8 we show that any bounded end of a rotating 
shrinking soliton must converge to these trivial cases, i.e. e eA invariant circles, or 
Abresch-Langer curves in N(A). 

The class of rotating shrinking solitons also includes the rotating solitons for 
Curve Shortening of immersed curves on the sphere § n_1 . These were studied by 
Hungerbuhler and Smoczyk in [10] (in [10] solitons on other surfaces were also 
considered). The connection is explained in §3.1 and §6.3. 

In Lemmas 8.1 and 8.3 we show that a shrinking-rotating soliton can have at 
most one unbounded end, unless it is a straight line in N(A). If a rotating-shrinking 
soliton has an unbounded end, then we find that there are two distinct cases. One 
possibility is that the distance from points on the end to N(A) never falls below a 
certain lower bound (specified in Lemma 7.7). In this case the end is asymptotic to 
a generalized logarithmic spiral. In Lemma 7.8 we prove existence of such solitons. 
The other possibility is that the unbounded end of the soliton stays within a fixed 
distance of N(A). In this case it must accumulate on Abresch-Langer curves and 
straight lines in N(A). Although we were initially inclined to conjecture that this 
last case could not occur, preliminary computations involving matched asymptotic 
expansions, which we intend to present in a forthcoming paper, indicate that that 
such exotic solitons may actually exist. 

Acknowledgments. The second author (SJA) began exploring the world of 
3D curve shortening solitons together with John Sullivan while at the IMA in 1990. 
This present paper owes much of its inspiration to those preliminary investigations. 

SJA was supported by NIH GM071794 and the Welch Foundation 1-1619. LFW 
was supported by NIH GM081549 and the Welch Foundation 1-1644. SBA was 
supported by NSF on grant DMS0705431. 

2. The symmetry group for Curve Shortening and its one-parameter 

subgroups 

The Curve Shortening equation is invariant under translations in space-time, 
rotations in space, and, being a heat equation, parabolic dilations of space-time. In 
this section we classify the one-parameter subgroups of the full symmetry group of 
Curve Shortening. We then describe what solutions invariant under each of these 
subgroups look like, and derive the differential equations they satisfy. 

Let & be the group of transformations of space-time M n x M of the form 

(5) g : R" x K -> R" x R, g(x, t) = (BDbc + a, 2 t + b) 

where 6 > 0, b € R, ft € 0(n,R), and a <E R n . Any transformation g e & leaves 
(1) invariant, i.e. it maps solutions of (1) to solutions of (1). 
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Let 

(6) 5e (x, t) = (0(e)K(e)x + a(e), 6(e) 2 t + b(e)) 

be a one-parameter subgroup of (5. Such a subgroup is completely determined by 
the vector field 

dg e (x,t) 



v<*t)= de 



= ((0 + M)x + v, 29t + w), 



e=0 

where 6,v,w, and M are defined by 
d_ 



(9(e), a(e), 6(e), 51(e)) = (9,v,w,M). 

o 

Specifically, the orbit <7 e (x ,io) = ( x ( e ))^( e )) through any point (x ,io) is the solu- 
tion of the first-order differential equation (x, i) — V(yi, t), where the dot denotes 
differentiation with respect to e. Thus, the general one parameter subgroup of is 
obtained by integrating the differential equations 

(7) x=(6> + M)x + v, i = 29t + w 

for suitable choice of constants 9 € R, w 6 R, v e R™ and M 6 so(n,]R) (i.e. 
M : R" -> R" is antisymmetric). 

Instead of considering the most general one-parameter subgroup generated by 

(7) , one can simplify the problem by conjugating the subgroup with some /i£0, 
i.e. by applying a linear change of coordinates (x, t) — h(x,i). Below we choose 
simple representatives in each conjugacy class of subgroups. 

Choice of parameter for the subgroup. If (9, v, w, M) generates the sub- 
group {g e }, then the one parameter subgroup {h e } generated by any multiple 
(X9, Av, Xw, AM) is given by h e = g\ e . Therefore multiplying the vector field V 
with a constant does not change the subgroup it generates. In particular, we may 
always assume that 

(8) w > 0. 

Rotating the coordinates. We can rotate in space and substitute x = Sx, t = t, 
for any orthogonal S e 0(n,M). Then t and t satisfy the same equation, while x 
satisfies 

dx/de = (9 + S^MS)* + S _1 v. 
Since M : R™ — > R™ is antisymmetric, its rank is even, say rank(M) = 2m, and 
there is an M-invariant orthogonal decomposition 

(9) R" = E 1 © E 2 © • • • © E m © N(M) 

in which each Ej is two dimensional. This allows us to choose the rotation S so 
that after conjugation with § the matrix of M has the form 

/ ui3 \ 



(10) M = 



with 3 = 



ijj m 3 

\ n _2m / 

where the rotation rates Wj are positive and ordered 

< w\ < uj 2 < ■ ■ ■ < 0J m , 
and where 0„_ 2 m is the (n — 2m) x (n — 2m) zero matrix. 



-1 

1 
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Cat. 


Parameters 




(*.*) 


ffe(x ,t ) 

= (x(c),t(e)) 


Description 


A 


6 = 
to = 
Mv = 


X 


= Mx + v 

i = 


x(e) = e eM Xo + ev 
i(e) = t 


Circles and 
Shrinking Helices 


B 


6 = 

MJ = 1 

Mv = 


X 


= Mx + v 

i = 1 


x(e) = e eM x + ev 
t(e) = t + e 


Translating and 
Rotating Solitons 


C 


= 1 
to = 
v = 


X 


= x + Mx 

i = It 


x(e) = e £ e eM x 
t(e) = e 2e t 


Dilating and 
Rotating Solitons 



Table 1. A taxonomy of one-parameter subgroups of 0. Up to 
coordinate changes there are three categories of subgroups of &. 



Translating in space. After rotating space, we can relocate the origin in space 
to p and substitute x = x + p, t = t. The equations (7) then become 

di/de = (8 + M)x + v, di/de = 28t + w 

where 

v = v + (8 + M)p 

If 8 ^ then, since M is antisymmetric, 8 + M, is invertible, and so we can choose 
p = — (# + M)~ 1 v to make v = 0. In short: if 8 ^ then, after choosing a different 
origin in space, we may assume that v = 0. 

If 8 — then v = v + Mp, so that we can still choose p to make v _L Range (M). 
Since M is antisymmetric, this is the same as choosing p so that M(v) = 0. 

Translating in time. Having chosen the origin in space, we can translate in time 
by setting t = $ + 1. This leaves the equation for x in (7) unchanged and turns the 
i-equation into 

di/de = 28i + w + 26d = 28t + to. 

When 8 = this has no effect at all, but if 8 ^ we can choose d = —w/28 so that 
w = 0. 

Dilating space and time. Finally, we can apply a parabolic dilation, x = t?x, 
t = i) 2 i, with d > 0. If 8 ^ then we have already arranged that v = 0, w = by 
appropriately choosing our origin in space-time. The equations (7) are in this case 
linear homogeneous and thus the dilation has no effect. 

If 8 = then the dilation turns the £-equation into dt/de = i)~ 2 w = w. Since 
w > by assumption (8) we are left with two options: either to > and we can 
arrange w = 1 by choosing d = to -1 / 2 , or to = 0. 

We conclude from the above discussion that after applying a suitable coordinate 
transformation we may always assume that the parameters 8, v, w, M belong to one 
of three categories, namely, A: 8 = w = 0, Mv = 0; B: 6 = 0, to = 1, Mv = 0; and 
C: 8 ^ 0, w = 0, v = 0. (See Table 1.) 
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3. Circles and Shrinking Helices 

Let {g e : e G R} be a one-parameter subgroup corresponding to category A and 
let {c(i) : t G R} be a family of curves invariant under g e . If (Co, to) is a point 
on one of these curves, then the g e orbit through (Co, to) lies in the curve. Since 
the g e do not change the time variable t, the g e orbit of any point sweeps out a 
curve at time t parametrized by e. Therefore a g e invariant family of curves can be 
parametrized by choosing one point Co(f) on the curve for each time t, and then 
letting the group g e act on that point. The result is 

(11) c(e,t) = e eM C (t) + ev. 

If the vector v/0 then one can always choose Co(t) _L v. 

Recall that we had chosen coordinates in which the matrix of M is given by (10). 
The exponential e eM therefore acts on R n by rotating the (x 2 j-i,x 2 j) component 
by an angle ujjt for j = 1, . . . , m. In particular, in two and three dimensions e e:M 
performs a rotation by an angle uiie around the origin in R 2 or the X3-axis in R 3 , 
and we see that c(e,t) traces out the circle in the Xicc2-plane through Co if v = 0, 
or a helix around the X3-axis if v ^ 0. 

In higher dimensions the curve traced out by c(e, i) depends on the rotation 
frequencies w\, . . . ,u> m . If they are integer multiples of a common frequency, i.e. if 
ujj = pjfl (f2 > 0, pj G N, 1 < j < m), then the exponential e eM is 2tt/H periodic. 
If v = then the curve c(e) = e eM C is a closed curve on the sphere with radius 
1 1 C 1 1 inIR 2m x{0} C R n . For instance, for n = 4, u x = pi, lj 2 = p 2 , with p lt p 2 G N, 
and for Co = (ri, 0, r 2 , 0) with r\ + r\ — 1, the curve 

c(e) = (n cospie, r\ sinpie, r 2 cosp 2 e, r 2 sinp 2 e) 

is a (pi,p 2 ) torus knot in S 3 . 

If the rotation rates u)j are not multiples of some common frequency, then the 
curve c(e) = e £M C is not closed, but instead is dense in a torus whose dimension 
is the largest number of ojj that are independent over the rationals. 

If v 7^ 0, then the projection of c(e) = e eM Co + ev onto R 2m x {0} is the curve 
c(e) = e e:vc Co. Therefore c(e) = e eM Co + ev describes a generalized helix over c. 

Suppose now that c(e, t) given by (11) is a family of curves evolving by Curve 
Shortening that is invariant under the one parameter group g e . Then Curve Short- 
ening implies a differential equation for the time evolution of C (t), which we now 
derive. 

Direct computation shows that 

c £ =Me eM C (t)+v, c ee =M 2 e eM C (i), c t = e tM C' (t), 

where ' is the derivative with respect to time. 

Since M is antisymmetric and since Mv = we have e cM v = v. From there we 
see that the two terms in c e are orthogonal. Therefore 

||c e || 2 = ||v|| 2 + ||MC (<)|| 2 . 
Hence the arc length derivative and derivative with respect to e are related by 
d _ 1 d _ 1 d 

ds ~ IM de ~ ^HvlP + IIMCoW*' 
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This allows one to compute 



eM M 2 C (t) 



||v|P + ||MCoWII : 



Together with c t = e eM C (i) this tells us that the family of curves (11) evolves by 
Curve Shortening if 

eM M 2 C (i) 



c t = c as , i.e. if e eM C (<) = 

ll v ll " 

Cancelling e eM on both sides we see that (1) implies 



HMCoWII 5 



(12) 



C'o(t) 



M 2 C (i) 



l|V || 2 + ||MC (i)|| 2 ' 

This equation is essentially linear, after a time-change. If one introduces a new 
time variable r, related to t by 

dt 



(13) 



dr 



\MC (t)\\ 



then Cq, as a function of r satisfies 



^° = M 2 C C 



Since M is given by (10), we have 



M = 

It follows that 
(14) 



-wile 



1 
1 



0/ 



c (r) = xy^Cj 



for certain constant vectors Cj G (£j as in (9)). To eliminate r, we observe 
that 



|MCo|| 2 = ||v|| 2 + 5> 2 e- 2 ^||C, 



and integrate (13) to get 
(15) 



t=\W\\ 2 r-J2 



— e~ 2 j t | 



One cannot explicitly solve this equation for t, but one can find approximate solu- 
tions when t — ¥ ±oo. Assume that lo\ < iqi < ■ • • < w m , and let fc and / be the 
smallest and larges values of j respectively for which Cj ^ 0. 

As r — > — oo we can approximate t by the largest term in (15), giving us 



t = -(- + <> 



-a-,r|| Cl ||2 



while, by (14), 



(1)) 
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These expansions for Co and t show that as r — > — oo, the actual time t also goes 
to — oo, and 



(i6) cm-/- 

Thus for t — > — oo the rescaled curves (— t) _1 / 2 c(e, t) converge to a circle with radius 
yj2/wi in the Ei plane. 

Examples: If C/ is the only non zero Cj, then (— i)~ 1//2 c actually coincides with 
this circle, but if there are other nonzero terms in (14), then when the Wj = pjfl 
are integer multiples of a common rate fi, the normalized curve (— t)" 1 ' 2 c(-, t) is 
a closed curve converging to a multiple cover of the circle with radius \f2uTi ; if at 
least two u)j with Cj ^ are independent over Q, then (— t)~ x / 2 c(-, t) is densely 
wound on a torus, which converges to the circle in Ei with radius ^/2uii as t — > — oo. 

The behaviour for t — > oo depends on whether or not v = 0. 

If v 7^ 0, then (15) implies that 

t= (||v|| 2 + (l))r, (r^oo). 

By (14), Co(t) decays exponentially as £ — !• oo, so that (11) tells us that the curve 
c(e, t) converges at an exponential rate to the straight line c(e, oo) = ev. 
On the other hand, if v = 0, then (15) implies that 



* = -(f +o(D) 
Reasoning as above we find that 



— Ul i T 

e k . 



< 17 » c »<^fS»§j + °(^'- 

Thus as t the rescaled curves (— i) _1 / 2 c(e, t) converge to a circle with radius 
y // 2/ujk in the Ef. plane. 

3.1. Connection with Curve Shortening on S" _1 . The solutions we have found 
for v = are related equivalent to certain special solutions of Curve Shortening on 

Lemma 3.1. If C (t) is a solution of (12) and ifv = 0, then ||C (t)|| 2 = 2(T-t) 
for some T € R. The corresponding curve c(t, e) = e eM Co(t) therefore lies on a 
sphere of radius \/2(T — t). 

Proof. We compute, using (12) and M* = — M, 

jn r m||2 _ o (C (t),M 2 C (t)) (MC (t),MC (t)) _ 

dt" ° W " ~ ||MC (t)|| 2 ~ IIMCoWH 2 ~ ' 

which immediately implies the Lemma. □ 

Lemma 3.2. Let c : (oo,T) x K — > R" &e a family of curves with \\c(t, e)|| 2 = 
2(T — f), and consider the family of curves c on i/ie unit sphere given by 



(18) c(t, e) = ^2(T - t) c(-i ln(T - t), e). 

Then c evolves by Curve Shortening if and only if c evolves Curve Shortening on 
sn-i 
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Proof. Since Curve Shortening for immersed curves in S™ is given by eg = c ss + c, 
the lemma follows from a direct computation: 

Ct = 7lN {c ^ e} ' Css = 7#=? 

where 9 = -\\n{T-t). □ 

Putting these two Lemmas together we see that if c(t, e) = e eM C (t) is a solution 
to (1), then for suitably chosen T the family of curves c(9, e) defined by (18) is an 
eternal solution of Curve Shortening on 8" . The asymptotics derived in (16) and 
(17) imply that c converges to a cover of some great circle on S™ _1 , both for 9 — > oo 
and for 9 — > — oo. 



4. SOLITONS 

The one-parameter symmetry groups that led us to the shrinking helices left the 
time coordinate t invariant. From here on we will consider one-parameter symmetry 
groups {g e } of © of categories B and C (see Table 1), for which dt/de ^ 0. 

In this section we will derive an ODE whose solutions generate all solutions 
to (1) that are invariant under a one-parameter subgroup of © of category B or 
C. In both cases we start with a family of curves c(£,i). Assuming the family is 
invariant under a category B or C subgroup {g e } we find an equation relating c t 
and c^. Combining this with the further assumption that c is a solution of (1), 
i.e. of p Ca c t = c ss , we can eliminate the time derivative c t after which we are left 
with an ODE. This ODE must be satisfied by the curve c(-,t) at all times t. 

4.1. The ODE for Translating and Rotating Solitons. If {g e } is of category 
B, then we have 6 = 0, w = 1, Mv = 0, and the subgroup is given by 

ff e (x,t) = (e eM x + ev,< + e). 

In words: g e advances time by e, and translates and rotates space by ev and e , 
respectively. 

If a family of curves c(£, t) is known at time t = t , and if the family is known to 
be invariant under the action of the subgroup g ei then this action determines the 
curves at other times. Thus for any given c(£o,£o) the point <7 e (c(£o,to),io) also 
belongs to the family of curves, so that a £i exists for which 

5e(c(£o,*o),*o) = (c(fi,*o + e),*o + e) 

holds. Different choices of e will lead to different values of £i. We denote the 
function £h>^ defined in this way by ^ = £(e). Consequently we have for all e 

(19) c(£(e), t + e) = e eM c(£ , t ) + ev. 

Since £(e) is obtained by solving this equation, and since £ <— > c(£, t) is an immersion, 
the Implicit Function Theorem implies that £(e) depends smoothly on e. We can 
therefore differentiate (19) with respect to e at e = 0, which leads to 

c ( + ^ct=Mc + v, 
de 

hence 

p Cs c t = p Cs (Mc + v) 
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and, because of (1), 

c ss = p Cs (Mc + v). 

This equation contains no time derivatives and is in fact an ODE. Therefore we 
find that at any time to the parametrized curve C(£) = c(£,£ ) is a solution of the 
ODE 

(20) C sa = p Ca (MC + v) 

Conversely, if a parametrized curve C : M. —> W 1 satisfies (20), then the family of 
curves defined by (19) satisfies (1). 

4.2. The ODE for Dilating and Rotating Solitons. If {g £ } is of category C 
then we have 9 = 1, w = 0, v = 0. The action of g e is given by 

5£ (x,;) = (eV M x, e 2e i). 

If c(£, t) is a solution of (1) that is invariant under {<? e }, and if it is defined at some 
time to € M, then, reasoning as in the category B case above, we find 

<7e(c(£ ,io),*o) = (c(£(e),e 2e i ),e 2e t ) 

for some £(e) that depends smoothly on e. Consequently, 

(21) c(e(e),e 2 %) = eV M c(£ ,t ). 

Once again we can differentiate with respect to e and set e = 0, this time with 
result 

2t c t + — c £ = c + Mc. 
de 

Taking the normal component, and using (1) we arrive at the ODE 

(22) 2t c ss =p Cs (c + Mc). 

We now have to distinguish among the possibilities to — and to =/= 0. 

If to = 0, then (22) is not a second order differential equation at all. However, 
(21) implies that for some smooth e H> £(e) one has 

c(£(e),0)=eV M c(£o,0). 

The curve at time to = is a reparamctrization of a generalized logarithmic spiral: 
it is traced out by e i— > e e e eM Co for some constant vector Co- 

If to 7^ 0, then we can divide both sides of (22) by 2to and use (1) to obtain an 
ODE for the curve at time t . Introducing the constants 

(23) a = and A = -^—M 

2to 2to 

we find that the parametrized curve C(£) = c(^,io) satisfies 

(24) C ss = p Cs (aC + AC) 

whenever c is a {g e }-invariant solution of (1). 

Conversely, if C : K — > E" is a given solution to (24) for some a^0 and some 
A E S0(M"), then upon defining to = l/2a, M = 2toA, one can use (21) to extend 
the curve C to a {<? e }-invariant family of curves c(£, t) that satisfies (1). We have 
not specified how to choose £(e) in (21), but in view of the invariance of Curve 
Shortening under reparametrization of the curves, the choice of £(e) is immaterial. 
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One could simply set £(e) = £o for all e, so that our extension through (21) of C 
would be given by 

c(£,t) = eV M C(0, with e= ^ log A 

Z to 

This family of curves is defined for all t that have the same sign as to- Thus if 
a given curve C satisfies (24) with a < 0, then the resulting family of curves is 
defined for all t < and represents a rotating and shrinking soliton. If a > then 
c is defined for alH > and represents a rotating and expanding soliton. 

4.3. The general soliton equation. Equations (20) and (24) are both special 
cases of (3), i.e. C ss = (a + A)C + v + AC S . Here A is a smooth function of 
the parameter and is uniquely determined by the requirement that C s _L C ss . 
Comparing the inner product with C s of both sides of (3) leads to 

(25) A = -((a + A)C + v,C s ). 
We can also write (3) in the following form 

(26) C ss = p T ((a + A)C + v), 

where T = C s is the unit tangent to the curve C. 

The soliton equation (3) defines a first order system of differential equations for 
(C,T), 

(27) C S = T, T s =p T [(a+A)C + v]. 

Since C G R™ can be arbitrary and T is always a unit vector, this system defines a 
flow on K" x § n_1 . It is this system that we used to produce the soliton images in 
this paper. 

5. Components in the Null space and range of A 

In this section we split K™ into the null space and range of A, i.e. R™ = N(A) © 
R(A), and we consider the corresponding components of a soliton. Thus, for a 
given soliton we write 

(28) C(s) = Z(s) + W(s), Z(s) e N(A), W(s) e R{A). 

It is well known that non-rotating dilating or non-rotating translating solitons 
(i.e. those with A = 0) are planar. This fact continues to be true in a weaker 
sense in the case where the rotation matrix A does not necessarily vanish. Here we 
show that the projection Z of any rotating-dilating or rotating-translating soliton 
onto the null space ^V(^l) is planar. In the case of purely rotating solitons, the 
projection Z is in fact at most a line segment, and the soliton is a graph over this 
line segment. 

Lemma 5.1. Let C be dilating-rotating soliton satisfying (3) with v = 0. Then 
the orthogonal projection of C onto the null space N(A) of A is contained in a two 
dimensional subspace of N{A). 

In particular, all purely dilating solitons (a ^ 0, A = 0, v = 0) are planar: they 
are copies of the Abresch-Langer curves. 
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Figure 2. Examples of non-planar solitons involving rota- 
tion. Shown are sixteen different solitons that rotate around the 
z-axis (top); rotate and translate upward along the z-axis (second 
row); rotate around the z-axis and expand (3rd row); rotate and 
shrink (bottom). 



Proof. C satisfies (3). Therefore, if IIjv denotes the orthogonal projection onto 
N(A), then Il^A — AUn — and thus Z = HnC satisfies 

Z ss = aZ + A(s)Z s . 

Let $i(s) and $2(5) be the two solutions of the general solution of the second order 
scalar linear differential equation $" = a$ + A(s)$' with $i(0) = 1, ^(O) = and 
$2(0) = 0, $' 2 (0) = 1. Then Z is given by 

Z(s) = $i(s)Z(0) + $ 2 (s)Z'(0). 

It follows that Z(s) is contained in the two dimensional subspace spanned by Z(0) 
and Z'(0). □ 

Lemma 5.2. Let C be translating-rotating soliton satisfying (3) with a = 0. Then 
the orthogonal projection of C onto the null space N(A) of A is contained in a two 
dimensional affine subspace of N(A). 
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In particular, all purely translating solitons (a — 0, A — 0, v ^ 0) are planar: 
they are copies of the "Grim Reaper. " 

Proof. We again consider the projection Z = n^rC. Since we always assume v £ 
N(A), we have ITvv = v. Therefore Z satisfies 

Z ss = v + A(s)Z s 

where s is still the arclength along the curve C. The general solution of the (scalar) 
ODE 

y" - \(s)y> = 1 

can be written as 

yOO=y(0) + i/ / (0)*iOO + * p (a), 

where $i is the solution of the homogeneous equation with <E>i(0) = and &[(0) = 
1, while $ p is the particular solution with $ p (0) = $p(0) = 0. 
It follows that 

(29) Z(s) = Z(0) + *i(a)Z'(0) + $ p (a)v. 

Therefore Z(s) is always contained in the plane through Z(0) spanned by the vectors 
Z'(0) and v. 

□ 

Corollary 5.3. Let T(s) be a smooth family of unit vectors, and let 

C(s) = C(0) + I T(s') ds' 
Jo 

be the corresponding curve in K. n . IfT(s) satisfies the differential equation 

(30) T s =p T (v) 

for some fixed vector v, i/ien the curve C is a purely translating soliton with velocity 
vector v. 

Lemma 5.4. Let C be a purely rotating soliton. Then its projection Z onto the 
null space N(A) is contained in some line £ C N(A). 

The projection consists of either one point, a half line obtained by splitting I at 
a point P € £, or the line segment between two points P,Q € I. 

If the projection is a half line or a line segment m C t, then the soliton is a 
graph over m in m x R(A). 

The projected soliton approaches any endpoint of m at an exponential rate, i.e. if 
linw ±00 Z(s) = P then ||Z(s) - P\\ < Ce~ 5 l s l for some C,S>0. 

This is illustrated in Figure 2 (top row). Non-planar rotating solitons in M 3 are 
graphs over the rotation axis, and if they are asymptotic to some plane orthogonal 
to the rotation axis, then they approach this plane at an exponential rate. 

Proof. By the previous Lemma the soliton is given by (29) where in our case v = 0. 
Thus the projection onto N(A) of a purely rotating soliton is given by 

Z(a) = Z(Q)+«i(«)Z / (0), 

where $i satisfies <&'/ = A(s)$' 1 with $i(0) = 0, $i(0) = 1. The projected soliton 
is therefore contained in the line I through Z(0) with direction Z'(0). If Z'(0) = 
then the projection reduces to just the one point Z(0). 
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Assume from here on that Z'(0) 7^ 0. The line I is parametrized by 

Z = Z(0) + 0Z'(O), <peR. 

The equation for $1 implies that $i(s) > for all s. One way to see this is to 
observe that $1 satisfies a linear homogeneous first order differential equation so 
that either <f>i(s) — for all s, or ^\(s) 7^ for all s. Alternatively, one can just 
solve the ODE for $1, resulting in 

(31) $i(s) = / S &fo' ' x(s " )ds " ds'. 

Jo 

Since $1 is a monotone function the map s i-> Z(s) is one-to-one. The range of $1 
is an interval (0_,<^ + ) where 

<j)± = lira $i(s) 

s— ^±00 

and the projection of the soliton is the line segment m C £ consisting of all points 
Z(0) + (^Z'(O) with </>_ < </> < </>_)_. For each € (0-, <fi+) one can find s = ^O) 
and write 

C(s) = Z(0) + ^Z'(0)+W, 

where W is the projection of C(s) onto TV (.A)- 1 - = i?(.A). Regarding W as a function 
of <f> £ (</>-, 0+) we see that the soliton is indeed a graph over the line segment m. 

We now estimate the rate at which $i(s) — > <j>± as s — > ±00 from (31), which 
requires us to look at A(s) given by (25). Under the assumptions of this lemma 
(25) reduces to A(s) = — (C S ,AC). Differentiate, use (3), and keep in mind that 
C s _L C ss : 

A'(s) = -(C SS ,AC) = -{C SS ,C SS - AC S ) = -||C SS || 2 . 

Thus A is strictly decreasing, except at points of zero curvature on the soliton. 

We distinguish between the following possibilities for the function A(s). 

A(s) = for all s. In this case A'(s) = and consequently C ss = for all s. 
The soliton must be a straight line. 

A(s) changes sign, i.e. A(s*) > > A(s* >t ) for certain s* < s**. Let A(s >tst ) = 
-S. Then A(s) < -6 < for all s > s,*. Hence $' 1 (s)/$' 1 (s„) = e^- A(s ' )ds ' < 
e -5(s-s»„)^ Integrating we find that <f) + is finite, and that |$i(s) — 4>+\ < Ce~ Ss for 
all s > s**. An analogous argument applies to the behaviour of 3>i(s) for s — > — 00. 

A(s) does not change sign, and A(s*) 7^ for some s* G M. Without loss 
of generality we assume that A(s) > for all s. As in the previous case where A(s) 
changed sign we find that $i(s) converges exponentially as s — > —00. We must still 
examine the behaviour of the soliton as s — > +00. 

The components Z and W of C satisfy linear differential equations 

(32) W ss -A(s)W s = .AW, 

(33) Z ss -A(s)Z s = 0. 

Multiply the second equation with the integrating factor and integrate. We get 

(34) Z s ( S )=e^o A ( s ') ds 'z s ( So ). 

We can choose so so that Z s (so) 7^ 0, for if we could not, then Z s (s) = for all 
s, implying that Z(s) is constant, in which case the soliton would lie in an affinc 
subspace parallel to R(A). 
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Since Z s is a projection of a unit vector, (34) implies that 

I X (s')ds' <-ln||Z s (s )|| 



for all s > so- The integrand A is nonnegative, so we can take the limit s oo to 

get 



(35) / A( S )d S <-ln||Z s ( So )||. 

J So 

Recalling that A(s) is non-increasing we conclude that lim^oo A(s) = 0. 
It follows from (34) and (35) that 

Z s (+oo) = lim Z s (s) 

s— ► oo 

exists, and is non-zero. 

Now that we know that A(s) = o(l) as s — > oo, we see that the equation satisfied 
by W is a small perturbation of the constant coefficient equation 

(36) W ss ~AW = 0. 

The eigenvalues of A, as a linear transformation on R{A) are {±mji, . . . ,±zw m }. 
It follows that the characteristic exponents of the equation (36) are 



±V±i^= (±l±i)y^p (fc = l, 



We note that none of these exponents lie on the imaginary axis so that the associated 
system is hyperbolic. A perturbation theorem [7] therefore implies that the general 
solution of (32) can be written as 

W(s) = W+(s) + W_(s) 

where W + (s) consists of exponentially growing terms, and W_(s) is exponentially 
decaying. 

We claim that the exponentially growing component W + vanishes. Indeed, 

A'( S ) = -||C SS || 2 = -||W SS || 2 -||Z SS || 2 

implies that f°° ||W ss || 2 ds < oo. Integrating twice, it follows that ||W(s)|| < Cs 3 ^ 2 . 
Since W and W + differ by the exponentially decaying term W_ we see that || W + 1 
also grows at most like 0(s 3 / 2 ). As W + grows exponentially this can only happen 
if W + vanishes. □ 

6. The Distance to the Origin and to the Null space of A 

6.1. A differential equation for distance functions. If C : K — > K™ is a soliton, 
i.e. a solution of (3), then we again split C = W+Z, as in (28). The lengths of these 
components satisfy various monotonicity properties, depending on the parameters 
a, A, v. 

Lemma 6.1. IfB : R" -> R n is linear with AT> = T>A, Sv = 0, then the quantity 

<S( S ) = i||03C|| 2 

satisfies 

(37) 6 ss -\(s)S s -2aS= pC^ 2 . 
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In particular, 5(s) = ^||W(s)|| 2 satisfies the differential equation 

(38) 8 SS - \(s)6 s - 2a6 = ||W S || 2 . 

If C is non-translating soliton, so that v = 0, then S(s) = |||C(s)|| 2 satisfies 

(39) 5 SS - \(s)5 s - 2a6 = 1. 

Proof. The general case implies the special case by choosing 23 to be the orthogonal 
projection onto R(A) for (38) or by letting 23 be the identity for (39). To prove the 
general case we simply compute the derivatives of 5 = |||23C|| 2 , 

S s = CBC.SC,), 
S ss = (23C S ,23C S ) + (23C,23C SS ) 

= ||23C s || 2 + (23C,23{(a + .A)C + v + AC s }) 
= ||SC s || 2 +a||SC|| 2 + A(SC,SC s ) 
= ||23C s || 2 + 2a<5 + Ac5 s . 

□ 

Lemma 6.2. Let C be a rotating soliton that is either translating (a — 0, Av = 0, 
v 0) or not shrinking (a >0, v = 0). Let 23 again be linear with AH — HA and 
23 v = 0. Assume moreover that i?(23) c R(A). 

If HC(s) ^ for some s G R, then ||23C(s)|| is either strictly monotone on the 
entire soliton, or else it has exactly one minimum. 

The same statement applies to the length o/W. 

In the next section we will use this lemma with 23 = A. In the case that has the 
most natural interpretation 23 is the orthogonal projection on R(A), so 23C = W. 

Proof. Since 23C(s) ^ for some s £ K, we have S(s) ^ for at least some seE, 
where 8{s) = i||23C(s)|| 2 . 

If 8 S 7^ for all s then S s must have the same sign for all s, and therefore 6 is 
monotone. 

The remaining case to consider is that S s {sq) = for some so- We will show 
that this implies S(s) has a strict local minimum at so, and in particular, that 6 
has no local maxima. This implies that a strict local minimum must be unique, 
since any two strict local minima would have to bracket a local maximum. Having 
established that S has exactly one critical point sg it follows that S(s) is decreasing 
for s < sq and increasing for s > sq. 

Let So be a critical point of 5. The differential equation for 5 implies that 

M*oH 2a<5(s ) + ||2C s ( So )|| 2 > 0. 

If either a5(so) > or 23C s (so) 7^ then we have S ss (so) > so that S has a local 
minimum at Sq. We are therefore left with the possibility that both ad(s Q ) ~ and 
23C s (s ) = 0. 

If a > then aS(so) = implies <5(s ) = and therefore 23C(s ) = 0. Applying 
23 to both sides of the differential equation (3) and using the assumption that 
23 commutes with A, we find that 23C satisfies a linear homogeneous differential 
equation, 



(40) 



23C SS = (a + yi)23C + A(s)SC„. 
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Figure 3. Plot of the distance to the rotation axis for a number of solitons. 

It then follows from BC(s ) = £C s (s ) = that BC(s) = for all s. We had 
assumed this was not the case. 

Finally we have one more case left, namely, a — 0. Compute the next two 
derivatives of 8 using (38): 

S sss = XS SS + X S S S + 2(23C S , HC SS ) 

S ssss — XS SSS + 2X S 5 SS + X SS S S + 2||SC S;S || 2 + 2(23C S , "BC SSS ). 

At so the conditions 23C S = and 6 S = 6 SS = imply first 5 SSS — 0, and then 
Sssss = 2||23C SS || 2 . 

If BC ss (s ) = 0, then (40) with a = implies that AT>C{s ) = SC SS - XT>C S = 
0. Since A is invertible on R(A) and since -ff(23) C R(A) this forces 23C(so) = 0. 
But then ¥>C{s ) = BC s (s ) = so that (40) implies T>C(s) = for all s, once 
more contradicting our assumption. 

Since we have shown that 25C ss (so) ^ Owe know that 6 SS ss(so) > and therefore 
S has a strict local minimum at sq is this last case too. 

□ 

6.2. Unbounded ends for non-shrinking solitons. 

Lemma 6.3. If C is a rotating non-shrinking soliton, i.e. a solution of (3) with 
v = and a > 0, then both ends of C are unbounded, i.e. 

lim ||C(s)|| = oo. 

s— >±oo 

Proof. We have shown in the previous Lemma that ||C(s)|| can have at most one 
minimum. This implies that either lim^oo |jC(s)|| = oo or else sup s>0 ||C(s)|| < oo. 
We will show that the second alternative leads to a contradiction. 

Suppose K — limsup^oQ ||C(s)|j < oo. Choose a sequence Sj — ► oo with 
l|C( s j)ll — > K, and consider the solitons defined by 
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On any bounded interval |s| < L the Cj are uniformly bounded. Their derivatives 
Cj iS are unit vectors and hence also uniformly bounded, and by virtue of the dif- 
ferential equation (26) the second derivatives are also uniformly bounded. After 
passing to a subsequence we may therefore assume that the Cj converge in C^ c to 
a solution C* of (26). This limit satisfies ||C*(s)|| < K for all s _ M, and 

||C*(0)|| = lim ||Cj(fl)|| = lim ||C( Si )|| = K. 

j— >OG j— >OD 

Thus C* is a rotating-expanding soliton on which S(s) = ^||C*(s)|| 2 attains a 
maximum at s — 0. By Lemma 6.1 we then have 

S ss (0) = \S S {0) + aS(0) + 1 = aS(0) + 1. 

Since a > we see that <5 SS (0) > 1 > so that 5(s) has a local minimum at s = 0. 
This is our contradiction. 

A slight modification of these arguments shows that the other end of the soliton 
is also unbounded, i.e. that lim.j_j._co ||C(_)|| — oo. □ 

Lemma 6.4. Let C be a translating-rotating soliton (a = 0. v ^ ; .Av = 0). Then 
z(s) = (v,C(s)) has at most one minimum, and 

lim z(s) = ±oo. 

s— ^_oo 

Furthermore, the R(A) component W (s) of C(s) satisfies 
(41) lim II W(_)|| = or oo. 

s— >±oo 



Proof. By Lemma 6.2 the quantity S(s) = ^||W(s)|| 2 has at most one minimum. 
It is therefore monotone for s sufficiently large. 
If 8 (s) remains bounded for s — > oo, then 

= lim 5(s) 

s— f oo 

exists. Consider the sequence of translates Cj(s) = C(j + s) — (j e N) 

of our given soliton. These curves are again solutions of (3) (since a = the 
equation is invariant under translations parallel to N(A)). Extracting a convergent 
subsequence one finds a soliton C*(s) for which 5l|W*(,s)|| 2 = Soo is constant. 
According to Lemma 6.2 this is impossible unless W*(s) = for all s. So we see 
that for our original soliton we either have ||W(s)|| — > oo or ||W(s)|| — > as s — » oo. 

The case s — > — oo follows by reversing the orientation of the curve and applying 
the same arguments. 

For z(s) one has a differential equation, 

(42) z ss = \{s)z s + ||v|| 2 . 

This equation implies, as before, that z is either strictly monotone or has a unique 
local minimum, which then also is a global minimum. 
We also have the following trivial bound for z s 

|2 S | = |<C S ,V)|<||V||. 

Case 1: assume that z s > for s > so The coefficient A is given by (25), 
which in the current context is A = (C S ,.AC + v). As before we have 

dX _ 
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We therefore know that A is decreasing. To show that z(s) — > oo we distinguish 
between a few cases, depending on the behaviour of A for large s. 

If A(s) > for all s > sq, then (42) implies z ss > \z s which, in turn, implies 

Zs(s ) J z s (s') J S0 
Since z s < ||v|| for all s, we find that 

r°° iivii 

(43) / X(s)ds < In -V^r < °o- 
The integrating factor for (42) is given by 

(44) m(s) = e- f °o Hs ' )ds '. 

In view of (43) we conclude that there exist m± such that 
< 77i_ < m(s) < m + for all s > So- 

Integrating (42) we get 

(toz s ) = m(s)||v|| 2 > m_ ||v|| 2 

m(s)z s (s) > to_||v|| 2 s — C 
=> z s (.s)>^||v|| 2 S -C. 

Thus z s — > oo, which cannot be, since z s < ||v||. The contradiction shows that 
A cannot stay positive for all s > so- Since A is strictly decreasing we must have 
A(s) < for all large enough s. We may assume that A(s) < — 5 is s > Si for certain 

Applying the variation of constants formula to (42) we find 

z s ( s ) = e /; i Hs ' )ds ' z s ( Sl ) + ||v|| 2 / e-^' ^ s "^ ds "ds' (drop the 1st term) 

Jsi 

> ||v|| 2 T e^' x ( s ") ds "ds' (use A(s") > A(s)) 

J si 

> l|v|| 2 [ S e^-^ds' 

1 _ P ^(s)(s-si) 

The formula (25) for A implies 

|A( S )| < \\AC(s)+y\\<C(l + s) 

since ||C( S )||<||C(0)||+| S |. 

Thus we have z s > C/(l + s), which, upon integration, shows that z(s) — > oo as 
s — > oo. We are done with the case in which z s > for s > s . 

Case 2: assume that z s < for s > so- If additionally we assume that 
A < for large enough s, then z ss = Xz s + ||v|| 2 > ||v|| 2 . Integrating, we find that 
z s would have to become positive for large enough s, which rules out the possibility 
that A(s) < for some large enough s. 

We are left with the case that A > for all s > so- Since A is decreasing this 
implies < A(s) < A(s ) for all s > s . 
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Again consider the integrating factor m(s) from (44). We have for sq < s < si 
(45) - m(s)z s (s) — -m(si)z s (s 1 ) + ||v|| 2 J m(s')ds' . 

J s 

All terms are positive so we can let s± — > oo and conclude that 

m(s)ds < oo. 



's 

Since m(s) is nonincreasing this implies that lim^oo m(s) = 0, so that we can also 
conclude from (45) that 

-z„(a) = ||v|l / -^rids'. 



m(s) 



Use < A(s) < A(sq) to estimate 



TO ( S ') = e ff\(s")ds" > e -A( 80 )(«'-s) 5 
m(s) 

and from there 

llvll 2 

-Z.{8) > 

for all s > sq- After integration this implies that lim^oo z(s) 



□ 



6.3. Rotating Shrinking Solitons in E™ and Rotating Solitons on S n_1 . We 
observed in §3.1 that there is a one-to-one correspondence between Curve Shorten- 
ing on the sphere S"^ 1 and certain shrinking solutions of Curve Shortening in R™. 
Under this correspondence some of the rotating shrinking solitons we study in this 
paper correspond to purely rotating solitons on S™ -1 as studied by Hungerbiihler 
and Smoczyk in [10] (e.g. see Figure 6 in [10]). In particular, every rotating soliton 
on S™ -1 is also a shrinking rotating soliton for Curve Shortening on W 1 . The follow- 
ing Lemma shows directly that there is a large class of rotating shrinking solitons 
that lie on a sphere. 

Lemma 6.5. Let a < and v = 0. Then the rotating soliton that is tangent to 
the sphere with radius (—a)" 1 / 2 is entirely contained in that sphere. 

Proof. The squared distance 5(s) — |||C(s)|| 2 to the origin satisfies (39). Unique- 
ness of solutions of differential equations implies that if S(sq) = l/(— 2a) and 
S'(s ) = holds for some s then 6(s) = l/(—2a) for all set. □ 

7. Behavior at infinity 

7.1. High and low curvature decomposition of soliton ends. In this section 
we consider the behaviour of solitons far away from the origin. The main finding 
is that when C (or rather its projection onto R(a + A)) is large, the soliton C 
decomposes into two parts. One part consists of short arcs of high curvature on 
which C is approximately constant. On such arcs the soliton is approximated by a 
solution of 

(46) C ss =p Cs (a) 

for some constant vector a (one has a « (a +A)C + v). The solutions to (46) with 
constant a are exactly the Grim Reapers with translation velocity vector a. 
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On the other part of the soliton the curvature C ss is small, and in fact much 
smaller than the remaining terms in (26). It follows that on these low curvature 
arcs the vector (a + A)C + v is nearly parallel to C s . Low curvature arcs are 
therefore approximately solutions of 

(47) C s = ±a(C) 
where 

(48) a(C) = (a + A)C + v, a=^-. 

Il a ll 

After reparametrization (47) is equivalent to a linear equation. Namely, if one 
introduces a parameter related to arc length via 

ds 



,4! ' " 1 Ha(C( S ))H' 

then (47) is equivalent with C a = ||a||C s = ±a(C), i.e. 

(50) ±^ = (a + A)C + v. 

da 

When a/0 we may assume that v = 0, and the solutions to (50) are given by 

(51) C = e ±ff ( Q+yl >Co 

for arbitrary Co- These curves rotate according to A and simultaneously move 
toward or away from the origin at an exponential rate e aa . In the simplest, two- 
dimensional case we get logarithmic spirals; in higher-dimensional cases we get 
curves that spiral inward or outward on cones. Below we show that the low curva- 
ture parts of any soliton with a/0 are actually 0(1) close to one of the generalized 
logarithmic spirals (51). 

If the dilation parameter a = 0, then we must allow v / 0, although we still 
have v G N(A). In this case we can split C into its components in the null space 
and range of A, 

C = W + Z, WeR{A), ZeN(A). 

Then (50) decouples into W CT = AW, Z CT = v, which shows that the general 
solution to (50) is 

(52) C = e^Wo + crv + Z 

with W G R{A), Z G N(A) arbitrary constants. 

The simplest case is again two dimensional with N(A) = {0} and v = 0. Here 
the rotating solitons are the Yin- Yang curves that spiral outward. On the other 
hand, the corresponding solution (52) to (50) parametrizes a circle; in particular, 
it remains bounded. This suggests that for a = the solution of (52) can pro- 
vide a good approximation of finite length segments of the solitons, but that the 
approximation is not uniform as a — > ±oo. 

7.2. Small and distant Grim Reapers. In what follows we will make the pre- 
ceding discussion precise. We write the soliton equation (26) as 

(53) T s = ||a|| p T (a), 
and consider the quantities 

M =(a,T), z, = ||a|| 4 (l-/?). 
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Here jj, the cosine of the angle between the unit tangent T and a(C). If /i = ±1 
then a = ±T. More generally, if /i is close to ±1 then a almost lines up with T. In 
particular, we have the identities 

(54) l|a-T|| 2 = 2(1- n), ||a + Tj| 2 = 2(1 + fi). 

Since C ss = a — (a, T)T = ||a|| (a — fiTj , and since C ss and T are orthogonal, one 
finds that the curvature of the soliton is given by 

(55) ||C„|| = ||a|| y/1 - /A 
Therefore 

v= (\\a\\\\C s . 



(56) ||a||-V s >l-M 2 - u-^sfi^- 



can be interpreted as a scale-invariant curvature of the soliton. 

We begin with an estimate for the rate at which fj, and v change along a soliton. 

Lemma 7.1. There is a constant C, which only depends on a and A, such that fi 
satisfies 

C 

||a|| "n s > 1-/T- J 

The quantity 

*, = ||a|| 4 (l-/x 2 ) 

satisfies 

— C 

(57) ||a|| _1 i/ s < — 2/ij/ + C\i\fv + - — —v when /i > 0, and 

\\a\\ 2 

— C 

(58) ||a|| _1 i/ s > — 2uv — C\i\fv — v when u < 0. 

Proof. We begin with 

a s = (a + A)T, a s = Ha]]- 1 p a (a s ) = \\al\- 1 p k ((a + A)T) . 
This leads to 

H. = (a,T s ) + (a s ,T) 

= ||a|| (a, p T a) + ||a||- x (t, p a [(a + A)T] ' 
= ||a||(l - /i 2 ) + Hair^PaT, (a + A)T) 

Next, we observe that ||p a T|| 2 = 1 — fi 2 so that 

||a||-V« > l-^-llair'Ha + ^llVl-M 2 

C 

= i- A1 2 -_yr^. 

The quantity v = ||a|| 4 (l — /i 2 ) satisfies 

||a||-V s = -2 M ||a|| 3 ^ + 4||a|| 2 ||a|| s (l - M 2 ) 

Using (56) and 

|||a|U<||a,|| = ||(a + yi)T||<||a + ^|| 
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we find that when \i > one has 

||a||-V s < -2 M ||a|| 4 (l- M 2 -^ v / r^) +4||a|| 2 || a + yi||(l 



_ C 
ll a ll 

as claimed. The inequality for [i < follows in the same way. □ 

The following lemma tells us where the regions of large curvature are. 

Lemma 7.2. There is a constant C < oo depending only on a and A such that for 
any point Cq — C(sq) on a solution of (26) that satisfies 

(59) A = ||a(C )|| > C, and |/i(C ,T )| < 1 - <?V 
there exist s_ < s + with sq G [s_,s + ] and 

(60) S+ -s_ < -r-lnA , 

such that at the ends of the arc C([s_,s + ]) one has 

(61) = -(1 - CA^), and - (l - CA^). 

Since the quantity fj, measures how well the unit tangent T matches the vector 
a(C), this lemma says that if at some point Co on a soliton the unit tangent is 
not very close to either a(C) or — a(C), then the point Co is contained in an arc 
on which the unit tangent T turns from almost — a(C) to +a(C). The larger A 
is, the smaller the upper bound (60) for the length of the arc; if the arc is very 
short then CwQ on the arc, and it approximately satisfies (46) with a = a(Co), 
i.e. the arc is approximately a Grim Reaper. See Figure 4. 

Proof. Consider the arc of length 2 centered at Co- Along this arc we have ||C — 
Co|| < 1, and hence 

ill - |K||| < ||a-a || < ||a+^l|| ||C-C || < C. 
If we assume that Ao > 2C then this implies that 

||a|| > A - C > \A and ||a|| < A Q + C < 2A . 
Along this arc it follows from (56) that 



AC 



2\ 



(62) ||a||-V s > 1 - A' 2 - ^42 v/W^ > |(1 - 

provided y/l — /x 2 > 8C/Ag, i.e. provided 

ImI < 1 - ca«\ 

Integrating the differential inequality (62) we find that the length of any arc con- 
taining Co on which the above inequality holds is at most 

l-CVt- 4 l-CA" 4 l-CA- 4 

f 2 dfi 4 f dn C, , 

" J ||a|| 1 - fx 2 - A J l-t^-A 

-1+CA- 4 -1+CA- 4 

If A is sufficiently large then CA^ 1 \nA < 1, and the arc is contained within the 
arc of length 2 with which we began the proof. □ 
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Figure 4. High curvature segments of solitons approximate Grim 
Reapers (Lemma 7.2.) The solitons here rotate around the z-axis 
with rate u>, and translate along the z-axis with velocity v z . Mag- 
nifications of the high curvature segments of the solitons are shown 
in the right column. 



Lemma 7.3. If ||a|| > C and n > 1 - C||a||" 4 , then 



(63) 



ds 



< CM 
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If ||a|| > C and /z < -1 + C||a||- 4 ,, then 
(64) rf||a| 



<C||a|| 



Proof. From a = [a + ^4)a + v we get 

^||a|| = (a,a s ) = (a, (a+A)T). 

If fi > 1 - C||a||- 4 , then (54) implies 

||a-T|| = v ^(r^)<q|a|r 2 , 

so that 



|a|| = (T,{a + A)T) + {k-T 1 (a + A)T) 



= a + (a - T, (a + A)T) 
where we have used that (T,AT) = and (T, T) = 1. Hence 
dlMI 



ds 



< |(a-T,(a + yi)T>| < ||a-T|| ||a + ^.|| ||T|| < C||a|| 



Thus (63) is proved. To prove (64) we assume fj, < — 1 + C||a|| 4 and use (54) again 
to conclude that ||a + T|| < C||a|| -2 . From there we get 

d||a|| 



ds 



< |(a + T,(a + yi)T)| < ||a + T|| ||a + yi|| ||T|| < C|ja|| 



□ 



Lemma 7.4. There is a constant C depending only on a and A such that for a > 
it follows from ||a|| > C that ||a|| is increasing when /i > 1 — C 1 1 a| | 4 and decreasing 
when fj, < — 1 + C||a|| -4 . 

On the other hand, if a < and ||a|| > C then ||a|| is decreasing when \i > 
1 — C||a|| _ and increasing when /i < —1 + C|ja|| . 

Proof. In the case a > we let C\ be the constant from (63). Then if ||a| > 
yJa/2C\ then (63) implies that <i||a||/<is > \a > 0. So the Lemma holds if we set 

The case a < is entirely analogous. □ 
Consider the two regions 

X+(K) = {(C,T 

= {( C ' T 

where a = (a + A)C. 

Lemma 7.5. For any C 6 R." one has 

\a\ ||C|| < ||(a + yi)C|| < \\a + A\\ ||C||. 

Proof. The second inequality is just the definition of the operator norm ||q + A\\. 
The first inequality follows from the fact that A is antisymmetric. By Cauchy's 
inequality and (C,yiC) = one has 

±a||C|| 2 - ±(C, (a + A)C) < ||C|| \\(a + A)C\\, 

which implies the Lemma. □ 



|a|| > K, v < K,[i > 0} 
|a|| > K, v < K, [i < 0} 
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fi=<a, T> 
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/i=<a, T> 






IHI n 


> 




I ll a 








/ 






1 L - 1 

a>0 


a 


< 



Figure 5. A caricature of the soliton flow (27) in the cases a > 
and a < 0. 

Lemma 7.6. //-ftT is large enough then for any solution (C(s),T(s)) to (^7j wii/i 
||a(0)|| > 2K one has 

• either (C(0),T(0)) £% + (K)uR-{K), 

• or there exist s_ < < s+ suc/i i/iai (C(s + ), T(s + )) G 3i + (_fC) and 
(C(s_), T(s_)) G 3?_(if). Moreover, s± are bounded by 

Isil^CKO^lnHaCO)!!. 
Proof. This is a reformulation of Lemma 7.2. □ 



Lemma 7.7. If a > i/ien for sufficiently large K the region 'R + (K) is forward 
invariant and the region %_{K) is backward invariant under the soliton flow (27) 
on K" x S™" 1 . 

// (C(s),T(s)) € ft+(C) /or aZZ s > ; then 

(65) C = e ,T ( a+ ^r + 0(e- a ' T ) 

for some constant vector T ^ 0, provided a is defined as in (49). 

If a < 0, then there is a K < oo such that any solution (C(s), T(s)) of (27) that 
remains in K_ (K) for all s > satisfies 

(66) C = e -"^+A) T + Q^aa) 

for some constant vector T 0, again, provided a is defined as in (49). 

Proof. In the region !R + (K) we have ||a|| > K, and hence, by Lemma 7.5, ||C|| > 
CK in R + (K). 

We also have v < K in R + (K). By definition of v this implies fi > 1 — Cif/||a|| 4 . 
According to Lemma 7.4 we find that for sufficiently large K one has 

(67) di|a|1 



a 

> - 
ds ~ 2 



in 3i + (K). Thus the only way that a solution could escape from H + (K) is by 
violating the inequality v < K . However, in 5l + (K) we have, by Lemma 7.1, 

, i dv 



ds 



< (-2/i + C , ||a|r 2 ) i / + C^ 



< 



(-2 + C||a|r 4 + C||a|r 2 ) I , + C^ 
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If at any point on a solution in H + (K) one has v — K, then one also has 
llair 1 ^ < (-2 + CK- 2 )K + CVK. 

If K is large enough this implies ^ < when v = K. Therefore "Jl+(K) is forward 
invariant when K is sufficiently large. 

The exact same arguments show that 5l-(K) is backward invariant, especially 
if one observes that the transformation (C, T) H> (C, — T) maps Oi + (K) to 3l_(K) 
and reverses the orientation of the soliton flow. Therefore forward invariance of 
R+(K) is equivalent with backward invariance of Dl_(JC). 

We now turn to the proof of (65). If a soliton satisfies (C(s),T(s)J € Oi + (K) 
for all s > 0, then we have v < K and hence /x > 1 — C||a||~ 4 for all s > 0. We are 
assuming that a > 0, so we have 

||a|| = ||(a + yi)C|| > o||C||, 

and therefore fx > 1 — C||C||~ 4 . Since ||a — T|| 2 = 2(1 — fj,) we also have 

(68) ||T-a|| < C||C||- 2 . 

Recall that a = a/||a|| and a = (a + A)C (since a > we have v = 0). We get 



\(a + A)C\\~ - {a+A)C 
as 



Defining a as in (49), we arrive at 
dC 



< ll(« + -A)C|| 1T ^ w < 



C 
ICI 



c 

icii 



(69) 



— -(a + A)C = g(a), where ||g(a)|| < ^ 



In particular, ||g(er)|| is bounded since ||C|| is bounded away from in 3l + (K). 
Multiply (69) with the integrating factor e a ( a + A ) and integrate: 



(70) 



e -^(«+^)c(^) = er^+^Cia) + f ' e- a ' {a+A) g{<7') da' 

J cr 



Since A is antisymmetric e aA is orthogonal and ||e cr '" + ' /l '' || — e aa . Therefore bound- 
edness of g(cr) implies that the integrand in (70) decays exponentially. This allows 
us to take the limit cr* — > oo and conclude that 



(71) T= lim e- a ^ a+A ^C{a,) 

cr* — f oo 

exists. Moreover, we have 

poo 

Y = e -^ a+A '>C(a)+ / e-°'( a+A ^{a')da' 

J a 

and therefore 



(72) 



C(a) 



_ a(a+A)-, 



g(a') da'. 



This last integral can be estimated by 



-(<r' -<r)(a+A) 



g(</) da' 



pOO 

< sup ||g(cr')ll / e- a ^'- a) da' 



- sup ||g(a')|| 
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In particular, the integral is bounded. It follows that r ^ 0, for if T were to vanish 
then (72) would imply that C(cr) is bounded for all a > 0. This is impossible when 
a > because (67) says that ||a|| must keep growing as long as the solution stays 
in 

Having established (72) with T ^ we can conclude that ||C|| > Ce aa . Com- 
bined with ||g|| < C/||C|| this implies ||g(<r)|| < CeT a ° ' . Our estimate for the 
integral in (72) then implies that this integral is 0(e _Q<T ), as claimed in (65). 

To prove (66) one can use very similar arguments. Since (C(s),T(s)) 6 Ol-(K) 
for all s > we have fi < -1 + C"||C||- 4 , and hence ||a + T|| < C||C||- 2 . From 
there one concludes that 

^ = -( a + yi)c + o(||c|r 1 ), 

da 

which then leads to (66). □ 

When a < the regions H±(K) are no longer invariant for the soliton flow in 
either forward or backward direction, and most orbits that enter D?_ (K) will leave 
31- (K) after a while. However, there do exist solutions of the soliton flow that, 
if you follow them in the direction of T, remain in 3i_(K) forever. We have just 
shown in Lemma 7.7 that such orbits are again asymptotic to logarithmic spirals. 
The following Lemma constructs these special solitons via a topological "shooting 
method." 

Lemma 7.8. Let a < 0. If ||a || > K then there is a T G § n_1 with ||T + 
a(C )|| < C^Haoll" 4 such that the solution (C(s),T(s)) of (27) with (C(0), T(0)) = 
(Co, To) remains in Jl-(K) for all s > 0. 

Proof. The map 4> : (C, T) H ► (a, T) = ((a + A)C, T) is a homeomorphism of the 
set D?_ (K) with 

- {(a,T) : ||a|| > K, (a, T> < -y/l - ^a]P}- 
The boundary of consists of two parts, 

K in = ^- 1 {(a,T) : ||a|| = K, (a,T) < -0 - X||a||- 4 } 

and 

K ex = ^.- 1 {(a,T) : ||a|| >JT, (a, T> = - ^1 - ^||a||- 4 } 
As above one concludes that within 3?_ (K) one has 

-f-||a|| > -a/2 > 0. 
as 

Therefore any solution that starts in 'R-(K) can never reach 3?; n again, since ||a|| = 
k on 3?j n ; if a solution starts in H_(K) and then exits 3?_(-£tT), it must hit 3? ox . 

The same arguments which used Lemma 7.1 to establish forward invariance of 
Jl+(K) when a > show in the case a < that, when a solution within 31- (K) 
reaches v = K, one has v s > 0. Hence solutions that start in and then exit 
Jl-(K) cross H ox transversally when they exit. It follows that for those initial 
data (C , T ) G R-(K) whose orbit {(C(s),T(s)) : s > 0} under the soliton 
flow eventually exits 3l-(K), the point -E(Co,To) G 3l ex at which they exit is a 
continuous (and even smooth) function of the initial point (Co, To). We call the 
map E : (C ,T ) H> E(C ,T ) the exit map. 

We now consider the topology of the exit set # ex . 
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For any fixed a with ||a|| > K the set of unit vectors T € S" 1 C R" that satisfy 

(a,T) = -Vl-^||a||- 4 

is an (n — 2)-dimensional sphere in a hyperplane in M." perpendicular to a. From 
this one sees that the map 

(?3) (C ' T) ^("'^) 

is a homotopy equivalence of 3i-(K) with the unit tangent bundle TiS" -1 of S n_1 . 
Under this homotopy equivalence the set of initial data (Co, To) with prescribed 
Co maps to the fiber over a . 

It is a fact from homotopy theory that the fiber in the unit tangent bundle 
TiS"- 1 S"" 1 is not contractible in TiS™" 1 . Indeed, the fiber bundle S"~ 2 
Ti§" _1 — » S™ _1 leads to a long exact sequence of homotopy groups 

■• • -► vr^^S"- 1 ) A vr n _ 2 (§"- 2 ) — ► 7 r n _ 2 (T 1 §"- 1 ) — ► 7r n _ 2 (S"- 1 ) -j- • • ■ 

Since 7r„_ 2 (§"- 1 ) is trivial tt^^S"- 1 ) is isomorphic to 7r n _ 2 (§™~ 2 )/d*7r n -i(§ n - 1 ) 
Both 7r n _ 2 (S n_2 ) and 7r„_ 1 (S™ _1 ) are isomorphic to Z, and by explicit computa- 
tion one finds that d* vanishes when n is odd, and acts by multiplication with 2 
when n is even. Thus for even n we find that 7r„_ 2 (TiS"^ 1 ) ~ Z 2 , while for odd n 
we have 7r„_ 2 (T'iS™ _1 ) ~ Z. In either case, the fiber S" -2 represents a homotopy 
class that generates 7r„_ 2 (S™ -2 ), and does not lie in the image of 9*. Its image in 
7r„_ 2 (T 1 §™ -1 ) is therefore nontrivial, and so the fiber is not contractible in TiEi 71 ^ 1 . 

To prove the Lemma we consider the given set of initial data: we are given one 
point Co, and we let T be any unit vector that satisfies 

(74) (a , T) < -v/l-^llaolh 4 

This set of points forms an (n — 1) -dimensional ball B' i_1 C %-(K). If the orbit 
starting at all possible initial data To in this ball B n_1 were to exit $.-(K), then 
the exit map E restricted to B™^ 1 would provide a continuous map B™ -1 — > 3i ex . 
On the other hand, the points on the boundary of B™ -1 are already in the exit 
set 3?ex, so the exit map restricted to <9B" _1 is the identity. Under the homotopy 
equivalence (73) the boundary <9B™ _1 maps to a fiber in the unit tangent bundle 
TiS" -2 . The exit map restricted to B™ -1 would then provide a contraction of this 
fiber within Ti§™ -2 to a point, which is impossible. We therefore conclude that for 
at least one choice of T satisfying (74) the orbit of the soliton flow will not leave 

□ 

8. Global behavior of Rotating-Dilating Solitons 

8.1. Decomposition of R™ into complex subspaces determined by A. In 

the previous section we have seen that rotating-dilating solitons can have ends that 
are asymptotic to logarithmic spirals. Here we determine all possible asymptotic 
behaviours that an end of a rotating-dilating soliton can have. 

There are a number of trivial rotating-dilating solitons, namely, circles with 
appropriate radius, and Abresch-Langer curves in case N(A) is nontrivial. We will 
show that all bounded ends of rotating-dilating solitons must be asymptotic to 
these simpler solitons. 
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It will be convenient to change our notation for the normal form of the infini- 
tesimal rotation matrix A. Recall that in a suitable basis A is of the form (10). In 
the most general case several of the rotation frequencies ujj may coincide. We now 
denote the set of distinct nonzero eigenvalues of iA by {±17^ : 1 < k < m} with 
< 17 1 < 0,2 < ■ ■ ■ < Ofn, and with fh < m. The matrix A can then be written as 



(75) A = 



( 17i2i 0\ 

OfnSrh 

V o ••• o i/ 



with 3 2 k =-l*v, 

and F k d = N(A 2 -Op). 



The subspaces F k are invariant under A. Restricted to each F k the matrix A 



satisfies A 2 = — 17|l, so that 



3* = nrU 



defines a complex structure on Fk- 

8.2. An almost Lyapunov function for the dilating-rotating soliton flow. 

The dilating-rotating soliton equation defines a flow on the unit tangent bundle 
1" x §" _1 of M™. When the rotation matrix A vanishes, one can use Huisken's 
monotonicity formula for Mean Curvature Flow [9, 5] to show that this flow is 
exactly the geodesic flow of the metric 

(d s ) 2 = e §H c ll 2 ||dC|| 2 

on W . For translating solitons one can use Ilmanen's monotonicity formula [11] 
instead to conclude that that self translating solitons are geodesies for the metric 

(d s ) 2 = e< v ' C >||dC|| 2 . 

Since geodesic flows are Hamiltonian systems one expects a high degree of recur- 
rence in the flow, and, indeed, all dilating, non-rotating solitons are either circles 
or Abresch-Langer curves. 

However, as we are about to show, when A ^ the flow becomes dissipative 
and most trajectories accumulate on a small number of circles or Abresch-Langer 
curves. Here we will exhibit an "almost Lyapunov function" for the flow (3) with 
A 7^ and v = 0. Recall that a Lyapunov function for a flow is a function which is 
non-decreasing along all orbits of the flow, and is in fact strictly increasing except at 
fixed points of the flow. A flow that has a Lyapunov function cannot have periodic 
orbits, and therefore the soliton equation cannot have a true Lyapunov function, 
since the dilating-rotating circles from §8.3 are periodic orbits for the flow. The 
following lemma presents us with a function that is non-decreasing along orbits of 
the soliton flow, and that is constant only on the trivial dilating-rotating solitons 
from §8.3. 

Lemma 8.1. If A ^ 0, then along any solution (C(s),T(s)) of the soliton flow 
the quantity 

V(C,T) - e* l|c||2+<v ' c) (T,yiC) 

satisfies 

(76) ^ = e fii c ii 2 +^ c )|| PT (yic)|| 2 . 

as 

In particular one has ^ > with equality only if AC is a multiple o/T. 
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Proof. Differentiating V, we get 



'~ = (aC + v,T)(T,.AC> + ~ 
as ds 

The last term can be expanded using antisymmetry of A and the equation for T s , 

d 

ds 



(77) e -fll c H -^ c ^ = (aC + v,T)(T,yiC) + ^(T,yiC) 

ds ds 



~(T,AC) = (T S ,AC) + (T,AT) 



e - f ||c||--(v,c>^l = ( aC + V;AC } + ||p T (yic)|| 2 . 



= (T S ,AC) 

= (p T (aC+AC + v),AC) 
= (p T (aC + v),yiC) +(p T (AC),AC) 
= ( P T(aC + v),^C) + ||p T (^C)|| 2 . 
Substituting this in (77) we get 

e -f ||c|| 2 -(v : c)^ = ( aC + v ,T)(T,yiC) + (pT(aC + v),yiC) + ||pT(^C)|| 
= ({«C + v, T)T + p T (aC + v),AC^ + ||p T (.AC)|| 2 . 

By definition of px one has w = Pt(w) + (w, T)T for any vector w. Therefore 

,dV 
ds 

The first term on the right vanishes. Indeed, antisymmetry of A implies that 
(C,AC) = 0, while Av = also implies (v,AC) = -(Av,C) = 0. This implies 
(76). □ 

8.3. The trivial rotating-dilating solitons. Under Curve Shortening any circle 
centered at the origin remains a circle with a varying radius. It therefore generates 
a dilating soliton, and if the circle happens to be invariant under the rotations e tA 
then it automatically defines a dilating-rotating soliton which satisfies (3). A circle 
can only remain invariant under the rotations e tA if it is contained in one of the 
eigenspaces F^. Moreover the circle must lie in a complex line in F/, relative to the 
complex structure 3k = A. 

If the null space of A contains a plane, then any Abresch-Langer curve in this 
plane defines a dilating soliton which also rotates according to A (by not rotating 
at all: e tA is the identity on the null space of A.) 

Lemma 8.2. Let a < 0, and let C be a rotating dilating soliton on which V(C, T) 
is constant. Then C is 

• a circle with radius l/y/—a in a complex line in one of the eigenspaces F^, 
or 

• an Abresch-Langer curve in the null space of A (or a straight line through 
the origin in N(A), which can be considered a "degenerate Abresch Langer 
curve. ") 

The functional V vanishes on the Abresch-Langer solitons. 
For the circles in Fk one has 

(78) V = ±-^, 

with the sign depending on their orientation. 
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Proof. Since V is constant on C we have ||yiC|| 2 = (T,AC) and hence AC = 
j(s)T for some smooth function 7. Keeping in mind that T is a unit vector, so 
that T, 1 T, we conclude from 



ds 



(AC, T> = (AT, T) +(AC, T s ) = 7 <T, T s ) = 



that 7 must in fact be a constant. 

If 7 = 0, then AC = 0. The soliton lies in the kernel of A, and must therefore 
be an Abresch-Langer curve. Clearly we have V = on any curve in N(A). 

If 7 7^ 0, then we have 



(79) 



d £= 1 -AC. 

ds 7 



On the other hand, by taking the inner product with T on both sides of the soliton 
equation (3) we get 

= (T,AC) + A(s)||T|| 2 = 7 + A(.s), 

where we have used C s = 7~ 1 j4C => C s _L C. It follows that A = —7. Using 
AC = 7T once more we can rewrite the soliton equation as 

C ss = (a + A)C + \T = aC. 

Differentiating (79) we also get 

C/g S — — — i/iC/g — — — — J\ CZ/. 

7 7 

Combining these last two equations we see 

A^C = 7 2 aC 

so that C must be an eigenvector of .A 2 with eigenvalue 7 2 a. Thus for some k we 
have C S Ffc and 7 2 a = — fi?, i.e. 



7 = ±" 



By equation (79) we have 



C(s) 



-as3, 



C(0), 



which shows that C is a circle in the complex line in that contains C(0). The 
radius of the circle is ||C(0)||, and is determined by the condition ||C S || = 1: 



ICJI = 



±V-a3 k e ± ^ s3 »C(0) = v/=5||C(0)|| = 1 

The value of V on this soliton is 

V^e-^ 2 (T,AC) - je- 1 ' 2 



|C(0)|| = 



1 



±- 



yj—ea 



□ 
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8.4. Compactiflcation of the Soliton Flow. One obstacle for studying the 
global behaviour of solitons is that the phase space R n x S ra_1 on which the soliton 
flow (27) is defined is not compact. To overcome this we show here that, after 
reparameterizing the flow, one can compactify the phase space and extend the flow 
by adding points at infinity. 

We map W 1 onto the interior of the unit ball via C H > P, 

P= y C . 

The inverse of this map is given by 

p 

C = 



v/1-HPll 2 ' 

A short computation shows that if (C,T) satisfy (27), then 

p a _ (c,c s )c 

S v/l + l|C|| 2 + 

This implies 



,3/2 • 



Vl-IIPPP, - (1 - ||P|| 2 )(T- (P,T)P). 
The evolution of T can be written as 

y/l - ||P|| 2 T s = p T ((a + A)y/l - ||P|| 2 C) = p T ((a + A)P). 

Hence, if we introduce a new parameter q along the soliton that is related to arc 
length via 

* = ; dS = VI + IICP da, 

then the soliton system (27) is equivalent with the following system for (P, T) 

(80) P, = (1-||P|| 2 )(T-(P,T)P) 

(81) T s = p T ((a + A)P). 

This system is obviously well-defined at all (P, T) g R™ x §™ _1 . The original phase 
space, which was described in terms of C and T, is embedded in the (P, T) phase 
space as {(P,T) : ||P|| < 1, ||T|| = 1}. The closure of this set, 

B" x S"- 1 = {(P, T) : ||P|| < 1, ||T|| = 1} 

is invariant under the extended flow (8.4), and it is also compact. 

The almost Lyapunov function V also extends to the compactified phase space, 
at least when a < 0. Indeed, direct substitution shows that 

V = e «ll p !! 2 /(i-ltP|l 2 )( 1 „ \\p\\2y 1/2 (AP,T) 

for ||P|| < 1. When a < the exponential is the dominant factor, which causes the 
whole expression to extend to a C°° function on B" x §™ _1 . One has V = at "the 
points at infinity," i.e. when ||P|| = 1. 

Since V is constant on the added points at infinity, one easily verifies that the 
continuous extension of V to B" x §" _1 still is an almost Lyapunov function for 
the flow (8.4): 




Figure 6. Rotating and shrinking solitons limit to a circle. 

Shown are solitons that rotate around the z-axis and shrink toward 
the origin simultaneously. A circle in the xy-p\a.ne with radius v2 
is one such soliton. Both ends of any generic soliton of this type 
are asymptotic to the special circle. 

Lemma 8.3. Let X(?) = (P(?), T(?)) be an orbit in B" x S™" 1 of the compactified 
soliton flow (8.4). Then q 1— > V{X{qf) is either strictly decreasing or constant along 
the entire orbit. 

If V(X(<r)) is constant then X corresponds to one of the trivial solitons from 
Lemma 8.2, or the entire orbit X(?) consists of "points at infinity," i.e. ||P(?)|| = 1 
for all cel. 

8.5. The ends of rotating-dilating solitons. Let C : K — > M. n be a solution of 
(3) with a < 0, and write X(s) = (P(s),T(s)) e B™ x S"" 1 for the corresponding 
solution to the system (8.4). By definition, the w-limit set of the orbit {X(s) : s e 
K.} consists of all possible limits X* = limn^oo X(s n ) one can obtain by choosing 
arbitrary sequences s n /* 00, i.e. 

w(X) d = (X» e B™ x S"" 1 : (3s„ /> 00) lim X, = X(s„)l 
= P| {X(,s') :s> >s}. 

The second description shows that w(X) is always a compact connected subset of 
B™ x §™ _1 , which is invariant under the flow (8.4). 
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FIGURE 7. A three dimensional projection of a shrinking rotating 
soliton in R 4 . Here a = —1, £l± = ^, fl 2 = 2. One end of the soliton 
is colored blue and the other red, so one can see that both ends 
are asymptotic to the same circle, but with opposite orientations. 
The limiting circle is the unit circle in the eigenspace F 2 (which 
appears as the larger ellipse in this projection), while the middle 
part of the soliton remains close to the unit circle in the eigenspace 
Fx. 

The a-limit set a(X) of the orbit {X(s) : s e R} is defined similarly as the set 
of all possible limits of X(s n ) for sequences s n \ — oo. 

The compactification of the soliton flow, and the almost Lyapunov function now 
give us the following description of the global behaviour of rotating-dilating solitons. 

Let C : K — > E™ be a rotating-dilating soliton, i.e. a solution of (3) with a < 0, 
and let X(^) = (P(<;), T(<j)) be the corresponding orbit of (8.4). Then, since V is 
monotone along the soliton, and since V is obviously bounded, the limits 

V ± = lim V(C(s),C s (s)) 

s— >±oo 

must exist. One has V- < V+. Moreover, V is constant on both w(X) and a(X) 
with 

V\ a( x)=V-, V\ U{X) = V + . 

By Lemma 8.3 both a and oj limit sets must consist of trivial rotating-dilating 
solitons, Abresch Langer curves in N(A), or points at infinity. 

We consider the possibilities for the w-limit set (the a-limits are analogous.) 

Lemma 8.4. If V + ^ then V + — ±fifc/\/— ea for some k, and w(X) consists of 
complex circles with radius (—a)" 1 / 2 in the subspace F^, as found in Lemma 8.2. 

In the special (but generic) case that all eigenvalues iClk of A are simple, the 
subspaces are two-dimensional, and they contain exactly one soliton- circle: when 
this happens w(X) contains that circle, and the soliton C is asymptotic to the circle. 
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The case V+ = is more complicated. If V+ = then w(X) could contain points 
at infinity, compact Abresch-Langer curves in N(A), or straight lines through the 
origin in N(A). 



Lemma 8.5. Assume V+ = 0. If w(X) is no£ contained in N(A) x S™ _1 then 
C(s) € 3^_(i4r) /or some large enough K and all large enough s. In this case the 
soliton C is asymptotic to a log spiral as in (66). In particular w(X) is disjoint 
from N(A). 

Proof. Since V+ = we know that w(X) is a subset of (<9B™ U N(A)) x S™" 1 . If 
there is a (Po,To) G w(X) n <9B™ with .APo 7^ 0, then we can choose a sequence 
s„ — > 00 for which ||C(s n )|| — > 00, and C(s n )/||C(s n )|| — > Po. 
In view of a(s n ) = (a +.A)C(s„) we have 

lim ||a(s n )|| = 00. 

If v(s n ) > K and if (C(s'J, T(s'J) g (if ) U (-*0 then Lemma 7.2 implies that 
for some = s„ + o(l) one has (C(s' n ), T(s' n )) G 3?+ (if). So after changing the s n 
slightly we may assume that (C(s„), T(s„)) G &_(if)u:R + (i<Q. If (C(s„), T(s„)) G 
3? + (if), then at s = s n we have 

(yic,T) = (yic,a) + (yic,T-a) 

_ [|(a+A)C|| h - 1 - a / 

' 4C||J 1 (>lC,T-a> 



|(a+A)C|| 



^^{^I^L-HT-all} 



On the other hand V(C(s),T(s)) is a nondecreasing function which converges to 
zero, so V(C(s), T(s)) < 0, and thus (AC,T) < for all s. Hence we have shown 
that at s = s n 

\\AC\\ IKo+^Cjj 

W - lien l|T ~ a|1 - l|a + yl " l|T_a|1 - 

Since we are assuming that (C(s n ), T(s„)) G 3?+(if) we know that ||T(s„) — 
a(s n )|| — > 0, and therefore we end up with 

C(.s„) 



AP 



lim v4t 



0. 



|C(s n ) 

This is impossible since we have assumed Po ^ N(A). The contradiction shows 
that we must have (C(s n ), T(s„)) £ for all n. Moreover, if (C(s„), T(s„)) G 

3?_(if) and if (C(s),T(s)) were to leave 3i_(if) for some s > s n then it follows 
from Lemma 7.2 that (C(s),T(s)) would end up in cR+(K). The previous argu- 
ments show that this again contradicts V+ < 0. We can therefore conclude that 
(C(s),T(s)) must remain in 3l_(K) for all s > s n . Lemma 7.7 implies that the 
soliton is asymptotic to a generalized logarithmic spiral. □ 

If w(X) is contained in N(A) x S™ -1 then all we can say now is that w(X) con- 
sists of lines and Abresch-Langer curves (and w(X) must of course be compact and 
connected.) It would be natural to conjecture that if w(X) contains an Abresch- 
Langer curve, then that curve is all of w(X), and the soliton converges to that 
Abresch-Langer curve. It is also natural to conjecture that w(X) can never con- 
tain a straight line, since the soliton would have to march up and down this line 
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infinitely often. However, preliminary computations involving matched asymptotic 
expansions seem to show that even in the 3 dimensional case a shrinking-rotating 
soliton exists for which a(X) is a circle in the xy-plane, while w(X) consists of the 
entire z-axis. The construction and analysis of this solution will be the subject of 
a future paper. 

References 

1. U. Abresch and J. Langer, The normalized curve shortening flow and homothetic solutions, 
J. Differential Gcom. 23 (1986), no. 2, 175-196. 

2. Steven J. Altschuler, Singularities of the curve shrinking flow for space curves, J. Differential 
Geom. 34 (1991), no. 2, 491-514. 

3. , Shortening space curves, Differential geometry: Riemannian geometry (Los Angeles, 

CA, 1990), Proc. Sympos. Pure Math., vol. 54, Amer. Math. Soc., Providence, RI, 1993, 
pp. 45-51. 

4. Steven J. Altschuler and Matthew A. Grayson, Shortening space curves and flow through 
singularities, J. Differential Gcom. 35 (1992), no. 2, 283-298. MR 1158337 (93b:58031) 

5. Sigurd B. Angcncnt, Shrinking doughnuts, Nonlinear diffusion equations and their equilibrium 
states, 3 (Grcgynog, 1989), Progr. Nonlinear Differential Equations Appl., vol. 7, Birkhauscr 
Boston, Boston, MA, 1992, pp. 21-38. MR 1167827 (93d:58032) 

6. Kai-Seng Chou and Xi-Ping Zhu, The curve shortening problem, Chapman & Hall/CRC, Boca 
Raton, FL, 2001. MR 1888641 (2003c:53088) 

7. Earl A. Coddington and Norman Lcvinson, Theory of ordinary differential equations, 
McGraw-Hill Book Company, Inc., New York- Toronto-London, 1955. MR 0069338 (16,1022b) 

8. Hoeskuldur P. Halldorsson, Self-similar solutions to the curve shortening flow, (2010), 
http://arxiv.org/abs/1007.1617vl. 

9. Gerhard Huiskcn, Asymptotic behavior for singularities of the mean curvature flow, J. Differ- 
ential Gcom. 31 (1990), no. 1, 285-299. MR 1030675 (90m:53016) 

10. N. Hungcrbiihlcr and K. Smoczyk, Soliton solutions for the mean curvature flow, Differential 
Integral Equations 13 (2000), no. 10-12, 1321-1345. MR 1787070 (2001j:53090) 

11. Tom Ilmanen, Elliptic regularization and partial regularity for motion by mean curvature, 
Mem. Amer. Math. Soc. 108 (1994), no. 520, x+90. MR 1196160 (95d:49060) 

DJA: St. Marks School of Texas, 10600 Preston Rd, Dallas, Texas 

SJA,LWu: Green Center for Systems Biology, UT Southwestern Medical Center, 
Dallas, Texas 



SBA: Mathematics Department, University of Wisconsin at Madison 



